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In this paper is discussed a class of static spherically symmetric solutions of the general relativistic 
elasticity equations. The main point of discussion is the comparison of two matter models given in 
terms of their stored energy functionals, i.e., the rule which gives the amount of energy stored in the 
system when it is deformed. Both functionals mimic (and for small deformations approximate) the 
classical Kirchhoff-St. Venant materials but differ in the strain variable used. We discuss the behavior 
". of the systems for large deformations. 
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I. INTRODUCTION 

Classical elasticity is a phenomenological theory to describe the properties of solids. It is used heavily in appli- 

cations such as structural mechanics and other engineering disciplines. This also explains the influence of classical 

Q elasticity in the formulation of the popular finite element method for the numerical solution of partial differential 
equations (see e.g. [1]). While it is entirely sufficient to restrict oneself to a classical theory for bulk matter to describe 
every-day engineering problems it is nevertheless of conceptual interest to formulate a theory of elasticity which is 
bJO compatible with the space-time structure established by Einstein's theories of relativity. 

' The first attempts to merge elasticity with special relativity go back to the early 20th century and there have been 

^ _ several other formulations including Synge [2] and Rayner [3]. The most influential work, however, has been the 
<^ . paper by Carter and Quintana [4] who formulated the geometric setting for the theory and derived the basic field 
equations. The theory has also been considered from a field theoretical point of view by Kijowski and Magli [5]. 
T— I Recently, the theory has been analyzed from the point of view of the initial value problem formulation by Beig and 
^-H Schmidt [6]. They showed that the field equations can be put into a first order symmetric hyperbolic form and they 
C5 prove among other things that the Cauchy problem for the system is well-posed under various circumstances. Based 
on this formulation it is shown in [7] that there exist solutions of the elasticity equations in Newtonian theory and in 
special relativity describing elastic bodies in rigid rotation. In [8] it is proved that there exist solutions of the static 
elastic equations for sufficiently weak gravitational interaction. Losert [9] analyzes the case of a self-gravitating 
elastic spherical shell and shows existence of solutions in the Newtonian case. 

In a series of papers jTollTTl [121 Karlovini, Samuelsson and Zarroug adopt the formulation of Carter and Quintana 
to discuss spherically symmetric equilibrium configurations and their radial perturbations. They also present an 
exact static and spherically symmetric solution with constant energy density. 

Our intention in this paper is to discuss two different equations of state in the static and spherically symmetric 
■ ■ context with two different materials. In [6] the familiar Kirchhoff-St. Venant stored energy functional for hyper- 
elastic isotropic materials has been extended to the relativistic case. Recall that this functional is quadratic in the 
strain variable and contains the Lame coefficients as two material constants. Kijowski and Magli in [5] use the same 
functional. However, they adopt a different definition for their strain variable which has a non-linear relationship 
to the strain used by Beig-Schmidt. Hence, this results in two different stored energy functionals which have the 
property that, by construction, they agree with the classical Kirchhoff-St.-Venant functional in the non-relativistic, 
small deformation limit. 

The plan of the paper is as follows. In sec. |ll] we provide the necessary backgroimd on the formulation of the 
theory of relativistic elasticity. The exposition follows that of L6J- We present the two formulations by Beig-Schmidt 
and Kijowski-Magli and point out their differences. It turns out that the only difference is in the definition of the 
strain variable which accounts for the above mentioned different energy functionals. 
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In sec.|lll]we specialize to the static and spherically symmetric case and derive the equations which govern this 
situation. We show that this system of equations has a unique smooth solution once the central compression of the 
body has been specified. 

Sec. |IV] is devoted to a study of various models. In order to compare the two energy functionals we consider 
various scenarios. We discuss a solid aluminum sphere and a relativistic highly compact material similar to the 
nucleonic matter inside a neutron star, described with both theories as well as with the classical theory of elasticity. 
We also discuss how the choice of different natural states affects the solutions. 



II. PRELIMINARIES 
A. Relativistic Elasticity 

The relativistic theory of elasticity in the form that we will use in this work has been described in 0] . The kinematic 
structure of the theory can be formulated as follows. As the basic variable one considers a (smooth) map 

f -.Ji ^ (2.1) 

from space-time {Ji, g) [23] to a 3-dimensional manifold the 'material manifold' or 'body manifold' or simply 'the 
body'. This is a reference manifold which carries some additional structure which will be described later. The body 
manifold can be interpreted as the collection of all point-like constituents (baryons) of the actual body. Coordinates 
on 3§ are labels for each individual 'particle' of the body. The map / is a map from a 4-dimensional to a 3-dimensional 
manifold so its derivative df must have a kernel. One requires that F = df has maximal rank at each point so that 
this kernel has dimension one and is spanned by a unit vector field n". Using small Latin indices for tensors on 
and capital Latin indices for tensors on the derivative of / may be written as = daf"^- In local coordinates {x") 
on ^ and (X"^ ) on ^ the map / is given by expressions of the form 

X^=f^{x"). (2.2) 

The dynamics of the theory is specified by a Lagrangian density p which is regarded as a functional of / and its 
first derivative df. In addition, it will depend on the metric g. Thus, the action may be written as 

<f,df;g] = J p[f,df;g] y^dit^d^x. (2.3) 

The Euler-Lagrange equations for this action are 

From the properties of / one can already derive several useful consequences. Let Cl^^c be a 3-form on This 3-form 
can be interpreted as defining a measure on ^ which gives to each subset of 3^ the number of particles contained in 
it. The pull-back of 0.abc to ^ along / is a 3-form tOahc on ^ which is dual to a vector field co". It is clear that this 
vector field spans the kernel of / so that it must be proportional to u". Hence we have the formulas 

^abc = ^ABcF^F^F^ = ^ahcd'^'^' 
, ,3 ^ ^abcd, , 

CO =nu = -e cvijcd- 
6 

The proportionality factor n is interpreted as the number density of particles (baryon density) constituting the body 
in the state it acquires when embedded into space-time. 

The energy-momentum tensor of the theory is defined as usual by the variation of the action with respect to the 
metric 

'^"'^ = + PSab- (2-6) 

A consequence of the diffeomorphism invariance of the Lagrangian is (see that 

V,T\ = ^aF^ (2.7) 
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i.e., that the elastic field equations are satisfied if and only if the energy-momentum tensor is divergence free. This is 
not necessarily the case in other field theories, such as e.g., for the Maxwell field. 

The inverse metric g"^ defines a contravariant, symmetric and positive definite 2-tensor H^^ on the body by push- 
forward with the map / (the minus sign is due to our signature) 

H^B ^ -g^bpApB^ (2.8) 

This characterizes the current state of the body which can vary due to the space-time curvature. In order to describe 
the variation the conventional way is to compare the actual state with a reference state that is given a priori as a 
fixed structure on the body This can be done by postulating the esxistence of a (positive definite) reference metric 
Jab on the body manifold which characterizes a 'natural' state of the body in which - by definition - there is no 
strain[24]. The difference E^^ = H"^^ — 7^^ between H^^ and the inverse 7^^ provides a measure of the 'size' of 
the strain on the body. Equivalently, one may use the linear map (f = E^g = ~ obtained by lowering an 
index on E^^ with Jab- 

Writing p = ne where e is the energy per particle then the second Piola-Kirchhoff stress tensor is obtained as the 
derivative 

de 

^AB=2^. (2.9) 

Thus, specifying e as a function of the strain provides the stress-strain relation i.e., the equation of state for the 
material under consideration. If the stress tensor tab does not vanish in the natural state in which there is no strain 
then one talks about a pre-stressed state, otherwise the state is called stress-free or relaxed. We will be concerned 
only with a relaxed state. Thus, the energy density has a minimum in the natural state. For most applications it is 
enough to assume that the energy density is at most quadratic in the strain and we will do so here. Invariance under 
coordinate transformations in the body implies that it can depend only on the scalar invariants of S" and since we are 
in three dimensions those invariants which are at most quadratic in S are Tr S" and Tr[S'^). Thus, the energy density 
can be written as 



p = nm[l+'^ |pTr(^2^ +,?(Tr^)2| 

\ + \ {piH^'^HAB - 2H^A + 3) + qiH^A - 3)'}) 



where m is the rest mass of a particle and p and q are constants. This is the stored energy functional which is 
assumed in [6]. It describes the so-called Kirchhoff-St. Venant materials. When we refer below to the Beig-Schmidt 
(BS) formulation we mean the use of this stored energy functional. 

The fact that there exists a metric on the body implies that there are now two 3-forms available: the 3-form CIabc 
which gives the number of particles in each sub-domain of the body and the volume form Vabc induced by Jab 
which gives the volume of the sub-domain. Since the two forms must be proportional we have 

Habc = noVABC, (2.10) 

thus defining the particle density ng in the natural state. This can be used to define the mass density po = '""o ir* the 
natural state. Using the 'natural' particle density Hq we can obtain the following formula 

^abc = ^ABC^a^b^c = noVABcF^Fj^ = nCahcdU'^ = n£ahc- (2.11) 

In local coordinates where CIabc = \/detj d^X and eabc = V— deth d^x with hah = Sab — Uaiib, we have 

Kov^deFydetF = nV-det/z. (2.12) 



B. The Kijowski-Magli strain 

The main difference between the Beig-Schmidt [6] and Kijowski-Magli [5] formulations is the choice of the variable 
which measures the deformation. Beig-Schmidt use the difference between the actual and the relaxed metrics on the 
body while Kijowski-Magli use a logarithmic variable. They claim that this variable has better behavior when large 
deformations are studied. 
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With our choice of conventions and notation this variable is 

S«'' = -^log(w«M^-7/), (2.13) 

where 7aj, = F^F;f 7ab is the pull-back of the reference metric on SS to the space-time. Note, that 7„;, is positive 
definite so that the tensor inside the parentheses has only positive eigenvalues and the logarithm is well-defined. 

Kijowski-Magli write down an action functional in terms of this variable. As before, the scalar character of the 
action implies that it can depend only on the scalar invariants of S and Kijowski-Magli assume that it is at most 
quadratic in S. They introduce the invariants 

a = TrS, ^^ilrS^ (2.14) 

where S is the trace-free part of S. Then, they write the action in the form 

^ = J n(^m + ^Ac(^ + B^^ ^y-detgd'^x. (2.15) 

Here, we have adapted the formula of Kijowski-Magli somewhat because we use the particle density n instead of the 
matter density and consequently we have to interpret e as the energy per particle. When we refer to the Kijowski- 
Magli (KM) formulation we mean the use of this stored energy functional. 

When deriving the equations of motion Kijowski and Magli use familiar techniques from Lagrangian field theory. 
However, their energy-momentum tensor is the canonical one and not the dynamical (symmetric) one which is ob- 
tained by varying the action with respect to the metric (see |13] for a thorough discussion of this difference). Since 
we are using the latter tensor we cannot simply take over the expression of Kijowski-Magli. Instead, we need to 
derive this energy-momentum tensor explicitly as given in appendix [A] We obtain 

Tab = P UaUb + noiAhab + nBSab- (2.16) 



C. Comparison of the two formulations 

In this section we want to compare the two presented formulations of relativistic elasticity. We establish that 
they agree on the linearized level and show how they differ for large deformations. In order to compare these two 
formulations we introduce the following variable which measures the deformation from a given state 

ej' = -{hac + lac)h'^ = -hJ'-^J'. 

In terms of ej' we can write the KM deformation tensor in the form 

S= -llog(l + e). 

The BS deformation is E^^ = H^^ - 7^^. We can relate these two difference deformation variables by the following 
computation 

Eab = E^-F^Eab = -h^'FfFilAClBDF^F^ - JacJbdJ^'' 

= -lacilbd + hdW^ = laceb- 

It follows that Ea^ = (H^^ - 7^^)7ab = eabh"^ = e/ and also Eab^^^ = eab^"^- Thus, the BS-energy density takes 
the form 



PBS = nmfl + ^ [peabe"^ + qie/f] 



The KM variables a. and /3 can be expressed in terms of e as well. Thus, e.g., a becomes 

1 

a = — - log det(l + e) 
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and, similarly, /3 can be expressed as before in terms of S and hence in terms of e. In order to connect with the BS 
formulation we expand the energy density up to quadratic terms in e. For the expansion of a and /5 we find 

so that the energy density of Kijowski-Magli up to second order in e is 

The expressions for the energy density in the two formulations agree in this approximation if we put 

1 

mp = B, mq — A — -B. 

The coefficients in front of the quadratic terms can be related to the classical elastic constants. Introducing the number 
density in the natural state one defines the Lame coefficients A = n^mq and Iji = nomp. Then, n^mA becomes the 
bulk modulus K while }i is the shear modulus of the material. 

Under these circumstances the energy-momentum tensor is given up to first order terms in eab by 

T^ab = —PO UaUb [Ijieab - Aec'Tafc) • 

We contrast this with the exact energy-momentum tensors 

n In 

BS: Tab = —po UaUb + -— {2}i 'yac£% + A ec^jab) (2.18) 
KM: r„fe = — pow«Wfc + — (2F/i«ba + AS„,,). (2.19) 
In the case of no deformation i.e., at a point where the body is in the natural state one has 

lab = -hah- 

The BS-energy-momentum tensor reflects this relationship. It almost agrees with the linearized energy-momentum 
tensor except that 7„;, appears instead of — Both theories are quadratic in their respective deformation variables 
and therefore describe in some sense a Hookean theory in which stress and strain are proportional. However, the 
relationship between the two different strain variables is highly non-linear. While the two energy-momentum tensors 
agree for small strain they disagree heavily for large deformations. Similarly, the stored energy functionals which 
give the energy per particle as a function of strain are completely different in the two cases when viewed in terms of 
the strain variable £a^- Thus, the two formulations describe materials with different equations of state. Both materials 
behave like the usual Kirchhoff-St. Venant materials for small strain, but have a completely different behavior for 
large deformations. We want to explore some of the consequences of these differences in the remainder of this article. 

In the KM formulation the strain variable is defined in terms of the difference tensors in space. This results in a 
term proportional to in the energy-momentum tensor, i.e. a term which is isotropic in space. In contrast, in [6] is 
used the difference between the actual and the reference state on the body as the basic variable. This results in a term 
proportional to 7„;„ i.e., isotropic on the body. This has the consequence that it is much easier to describe a fluid as a 
special case of elastic material within the KM framework than in the BS case. 

III. SPHERICAL SYMMETRY 

Now we specialize to spherical symmetry. We take the space-time metric to be the general spherically symmetric 
and static metric 

g = £^'7 dt^ - e^? dr^ - ? {de^ + sin^ 6 d<p^) (3.1) 
and we assume that the body metric 7 is spherically symmetric as well] 25], i.e., when expressed in polar coordinates 

7^B = e2?o rfj^2 _^ ^2 i^^qI ^ gj^2 ^$2) _ (3 2) 
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Since the geometry at the origin should be regular we need e^^" = 1 at the origin, i.e., fo(0) = 0. The function f and 
t] depend on r while depends only on R. The map / : ^ — > ^ is assumed to be equivariant and thus without loss 
of generality it can be expressed as 

f{r,6,<p) = {F{r),e,<p) (3.3) 
for some function F(r) with F(0) = 0. Then, the deformation gradient is given as 

= F'dr®dR+de®dQ + d(p®d^. (3.4) 
Clearly, because of staticity we must have m" = e~''3f and from (|2.12l l 



^=F'^^^ = F'.fo-4. (3.5) 
no etr^smfl 

The pull-back ^ah of the reference metric on 3S is 

7„,, = {e^°¥'fdr^ + {dd^ + sin^ 6 d(p^) (3.6) 
With these formulas and the abbreviations x = F'e^f^^^"' and y = F fr we can compute the deformation tensor 

= (x^ -l\dr®dr+ (y^ - l) {dO (E}dg + dcp(E} d^) (3.7) 



Using this variable and the formulas (|2.181 and (|2.19l l we can find the energy-momentum tensors in both cases. They 
are given explicitly rn the appendix |Bl 

A. The equations 

The Einstein equations in the spherically symmetric and static case are well known, see e.g. lO . They are 



Go° = e-'^[^ (l - e^^) - = -871 (3.8) 
Gji = (l - e^i) +^A= 871 P, (3.9) 



r 



= \n" + {rif - r^l! n } = 871 Q, 



(3.10) 



where we have put Tq^ = p, P = — Tj^ and Q = — T2^. These are three equations for the unknown functions rj 
and / (the function which specifies the reference metric is considered as given). A consequence of the Einstein 
equations is that the divergence of the energy-momentum tensor vanishes identically 

VaT\ = 0. (3.11) 

Under the current conditions this equation has only one non-trivial component 

P' + t]'{p + P) + -{P-Q) = 0. (3.12) 
r 

In order to obtain a useful system one replaces the equation (|3.101 by l |3.12|l . Furthermore, one integrates (|3.81 by 
introducing the mass function 



M(r) = An / pr^dr (3.13) 
Jo 

or, equivalently, the mean density lu = M{r)/r^ to obtain 

e-^^ = l-2r^iv{r). (3.14) 
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Inserting this into | |3.9l l one can solve for r]' and insert this into ( |3.12| |. Then the following system of equations is 
obtained 

rw' = -3w + Anp, (3.15) 

, 2^71? + w 

''1 = 1 — 3.16 
1 — zr-^zf 

rP' = -r^^^^ip + P)- 2(P - Q). (3.17) 
1 — zr^w 

This system is somewhat deceptive, because p, P and Q are functions of F and its derivatives. Since they contain F 
and F' in a non-linear way the third equation gives a complicated non-linear equation for F". Equivalently, we will 
regard these functions as depending on x and y defined above. Then P' — Pxx' + Pyy' ■ From their definition we get 
a relationship between x and y 

ry' = xe?-?o - y (3.18) 
which can be used to substitute for y'. With this preparation we now have the following final system of equations 

rw' = —3w + inp, 

ry' = — p~^° — y, 

^ Vl - 2r2a; ^ (3.19) 

j_Py(,, ^ „-fo^ ^^'^nP + w p + P 2 



Once a solution of this system is found we can obtain rj by integrating | |3.16l l, e~^^ is given by | |3.14| | and F is found 
from the definition of x. The functions p, P and Q are specified by the choice of the elastic model as functions of x 
and y, while e^^o ig any given function of r characterizing the natural state of the body. It is only restricted by having 
the value of unity at the origin. 



B. Behavior at the origin 

The equations are singular at the origin r = and it is not a priori clear whether there exist regular solutions. If 
there are solutions which are bounded near the origin then they have specific values there which can be obtained 
from the system by putting r = 0. Then the left hand sides vanish and from the right hand sides we get 

iv{0) = ^p{0), y(0) = x(0) =: fl, P(0) = Q(0). (3.20) 

This shows us that the only free datum is the value a. It characterizes the volume change of the body at its cen- 
ter. Since the body should be compressed we assume that a = lim,-^of ('')/'' > 1- The initial value for iv can be 
computed from the expression of p in terms of x and y. The third condition states that in the center the radial and 
the tangential stresses should be equal and this is a condition on the matter model which cannot be influenced by 
specifying initial conditions. The fact that the central compression is enough to characterize a solution uniquely is 
physically reasonable and corresponds to the fact that a static fluid configuration is uniquely characterized by the 
central pressure. 

In order to show that with these initial conditions there exist regular unique solutions we apply the theorem by 
Rendall and Schmidt [15]. The verification of the conditions necessary for that theorem are somewhat lengthy and 
we refer the reader to appendix [C] It follows from this analysis that for a given value a there exists a unique and 
smooth solution of the system of equations | |3.19l l in a neighborhood of the origin. 



IV. NUMERICAL MODELING OF SPHERICAL ELASTIC BODIES 



A. The models 



In the rest of this paper we solve the system I I3.19I I for several specific matter models. We consider two situations, 
a sphere consisting entirely of an ordinary material such as aluminum and a sphere which consists of material which 
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resembles the neutron star crust. In both cases we choose the two different energy functionals corresponding to 
the BS and KM formulation, respectively. For aluminum we use the values pQ = 2720 kg/m'^, A = 63.3 GPa and 
ji = 25 GPa. 

For the neutron star matter we follow the presentation in |16] where the structure of the neutron star crust is 
described in detail. In the crust of a neutron star the density increases from the outer layer with lO^g/ cm'^ to the 
inner edge where the density is approximately lO^^g/ cm'^. While the ground state of the matter is a lattice which has 
anisotropic elastic properties it is customary to approximate it by a homogeneous and isotropic elastic material. This 
material is under high pressure and hence it is much easier to shear it than to compress itl,26]. In fact, the material is 
often assumed to be incompressible. The shear modulus of the matter in the neutron star crust has been calculated 
in e.g. [17] and we use a value of = lO^^Pa. The fact that the crust material is almost incompressible means that 
the bulk modulus K is very large compared to the shear modulus and hence that also the Lame coefficient A is very 
large. We take it here three orders of magnitude larger than the shear modulus, i.e., A = lO^^Pa. 

In the relativistic theories there is no canonical choice for the relaxed state of the body. While it seems natural to 
specify a flat metric on ^ this is not necessary. The choice of the relaxed metric has been discussed in the litera- 
ture OH^. We follow here a suggestion by Carter li20ll according to which one can obtain the relaxed state by the 
following procedure. One assumes the body is heated up until it melts and then one lets it cool down imtil it solidi- 
fies again. Assuming that the fluid phase is an ideal fluid then the body settles in a state which can be described by a 
solution of the perfect fluid equations. Hence, besides a flat metric we also consider the spatial metric corresponding 
to an incompressible fluid with a constant density po, i.e, we put 

.-2?o = 1 _ ^^p,r\ (4.1) 

As a third formulation we consider the classical non-relativistic theory of elasticity. The equations for the classical 
theory can be obtained from the relativistic equations as the Newtonian limit, see |6]. The difference to the relativistic 
equations is that one puts ^ = = 0/ so that x = F'. Furthermore, eq. ( |3.12| | is replaced by 

P' + ri'p + ^{P-Q)=0, (4.2) 
where tj' is the gravitational force, determined from the equation 

V' = ^= rw. (4.3) 



The stress components P and Q have the same functional form in terms of x and y as those for the BS-energy- 



momentum tensor, while p is the mass density in the actual state, given hy p = jj-pQ = xy^pQ. 



Hence, the non-relativistic system is 

riv' = —3w + inp, 
ry' = x- y, 



rx' = ^{y-x)-'^wp-hp-Q). 



(4.4) 



All the numerical solutions have been obtained using the Runge-Kutta ODE solver suite provided in MATLAB. 
The calculation is started with an initial value a for x(0) = i/(0) which is used to calculate the initial value for 
iw(0) = 47r/3jO(0) from the energy-momentum tensor. The calculation stops when P vanishes, indicating that the 
boimdary of the body has been found. 



B. Numerical examples 

We first study the aluminum sphere for the three formulations of elasticity. Clearly, for small values of the relative 
central compression 6 — a — 1 the three formulations should be almost identical. We show in Fig. [1] the profile of 
the average density across the sphere for the three formulations for 5 = 0.001. In this case we find a sphere with a 
radius of 535 km and a mass of 1.74 x 10^^ kg. The figure shows the relative difference Aif = {w — wj^i) / Wf^ between 
the classical solution w^q and the BS and KM solutions w, respectively. The BS solution is indistinguishable from 
the classical solution, the maximum value of the relative difference being 6 x 10^^^, while the KM solution already 
indicates its general property: the system is more tightly bound than in the classical or BS case. Still, in this situation 
of small relative central compression the maximal difference is only 7.5 x 10~^. 
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FIG. 1: Relative difference between BS resp. KM and the classical solution for an aluminum sphere with relative central compres- 
sion of <5 = 0.001 



1. Aluminum with BS -formulation 

Let us now look at the BS model in more detail. The radial pressure is given by (|B2|| 

P = i(xy2)x2 ((A + 2^) (x2 - 1) + 2A(y2 - 1)) . (4.5) 

On the boundary of the body, this expression vanishes. This can happen either when x = or y = or if the term in 
parentheses vanishes. However, it follows from the equation for y' that as long as x remains positive we have y' > 
for y = so that y cannot vanish before x vanishes. Thus, on the boundary we have either x = or the pair (x, y) 
lies on the ellipse defined by 

(A + 2^)x^ + 2Ay2 = 3A + 2pi. (4.6) 

In Fig.|2]we show a sequence of such final pairs (x, y ) obtained from initial values S in the interval [0.01, 4] . Obviously, 
both cases discussed above can occur. For small relative central compressions the final pair (x, y) lies on the ellipse 
and for increasing compression it moves towards the y-axis until it hits it for an initial value of (5 ~ 0.585. Then it 
moves along the y-axis for unlimited values of y. The vanishing of x at the boundary means that the radial distance 
between two adjacent particles there becomes infinite, i.e. the body ruptures. Imagine a large elastic sphere without 
gravitational self-interaction being compressed so that the central compression is above the critical value. When 
gravity is switched on, the sphere will be divided into a central piece and a shell at the radius where x vanishes. 

The equation for x' in ll3.19l l shows that x vanishes with an infinite negative slope because the leading term on the 
right hand side goes like 1/x near x = 0. Thus, the solution becomes singular just at the boundary. 

The two different cases just discussed can also be seen in the behavior of the mass-radius diagram in Fig.|3l where 
we display radius and mass of the aluminum spheres corresponding to relative central compressions S G [10"'', 10'^]. 
We plot it in double logarithmic and linear axes. The curve shows three different regimes, the classical one where 
M <x (indicated by the solid line) and an 'extreme' regime where M oc R^, indicated by the dashed line and 
finally a 'linear' regime with M (x R where mass and radius decrease with increasing central compression. The 
cross indicates the configuration which is closest to the critical configuration where the radial strain x vanishes. 
This mass-radius diagram should be compared with Figure 1 from [12]. The similarity of the qualitative behavior is 
obvious. Karlovini and Samuelsson argue that the branch from the maximal mass towards zero is unstable and we 
do find numerical indications of this here as well. Increasing the central compression beyond the value needed for 
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FIG. 2: A sequence of pairs {x, y) for relative central compressions 0.01 < J < 4 
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FIG. 3: Double logarithmic (left) and linear (right) plot of the mass-radius diagram for BS aluminum spheres with relative central 
compression 5, logjg & € [—3, 3], in units of the earth mass and radius, respectively. The solid line is M (x J?^ 



the maximal mass configuration we observe that we can generate the smaller configurations up to a certain value 
of 5 depending on the required precision. Beyond this value the solver suddenly settles to a solution which yields 
a configuration in the 'eye' inside the mass-radius diagram. This dot in fact contains nine different configurations. 
The location of the 'eye' is roughly at the mass resp. radius for which the radius resp. the mass are maximal on the 
curve. The behavior of this system close to the eye should be analyzed in much more detail using more accurate 
solution methods. 



2. Aluminum with KM formulation 



In the formulation of Kijowski-Magli the radial pressure P is given by dBSI l 

P= (xy2)((A + 2;i)logx + 2Alogy). 



(4.7) 
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As before, at the boundary we have either x = or 

(A + 2/i)logx + 2Alogy = 0, (4.8) 
the case y = being excluded as before. However, now, a final pair needs to lie on the curve defined by 

1 2u 

This curve approaches the y-axis but never intersects it. This indicates that only the case when the final pair lies on 
the curve does occur. This is in fact confirmed in Fig. |4] where we show the final pairs (x, y) for aluminum spheres 

8 
7 
6 
5 

>,4 
3 
2 
1 



-0.2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

X 

FIG. 4: A sequence of pairs {x, y) for relative central compressions 0.001 < 5 <25 

with relative central compression 5 G [0.001,25]. While in the BS case the value of y can grow arbitrarily, this is not 
the case here. In fact, the numerical investigations show that the exhibited value of y w 4.1 is the maximal value that 
y can achieve. This behavior can be imderstood when we show the mass-radius diagram for KM aluminum spheres 
in Fig. |5]which shows a peculiar spiral. The maximal value of y is reached at the same point as the maximal radius. 
Thus, it is not possible with the KM formulation to create arbitrarily large objects. There exists a maximal mass and 
a maximal radius for KM aluminum spheres achieved for different objects and there exists a region where a KM 
aluminum sphere of a given radius can have at least four different masses. It looks like the sequence converges to a 
limit point. We have not been able to prove this rigorously. 

3. The neutron star matter 

We have also looked at an exotic material which is somewhat similar to the nucleonic matter that is assumed to be 
present in neutron stars. We show in Fig. [6] and in Fig. [7| the mass-radius diagrams for the neutron star like matter 
distributions with the BS and KM stored energy functionals. In both cases the diagrams look qualitatively the same 
as those for aluminum except that the size of the configurations are orders of magnitudes different. In the KM case 
we find a spiral as before while in the BS case we have the 'loop' with a linearly decreasing branch. Again, this 
branch seems to be imstable and the final dot in the diagram is the last for which we could generate a configuration. 

This shows that there is no qualitative difference in the behavior of aluminum and the exotic matter. This might 
change if one would use the high-pressure formulation developed by Carter and Qurntana |^]. 
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FIG. 5: Mass-radius diagram for KM aluminum spheres with relative central compression 5, logjQ 5 6 [—3,3], in units of the earth 
mass and radius, respectively. The solid line indicates the curve M rx R^. 




FIG. 6: Mass-Radius diagram for BS neutron star like matter with relative central compression 5, log^^Q 5 6 [—3,0.114], in units of 
solar mass and kilometers. 



4. The role of the relaxed metric 

As discussed above we employ two possible choices for the metric of the relaxed state of an elastic configuration. 
To compare the two different scenarios we compute configurations with the same relative central compression 5 for 
values of 6 between 10^^ and 1 for the two energy functionals and the two possible materials. In Fig.|8]we show the 
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FIG. 7: Mass-Radius diagram for KM neutron star like matter with relative central compression S, logjQ S 6 [—3,3], in units of 
solar mass and kilometers. 
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FIG. 8: Maximal relative difference in the mean density w between calculation with flat and curved metric for the relaxed state 
for given S with logjg S 6 [—3,0]. 



behavior of the maximal absolute value 

Aw 

= max 

IV r 

of the relative difference of the mean densities ivj- and Wc for the flat and curved cases resp. as a function of S. 
Obviously, in the given range of S the difference between the two configurations is almost negligible. The difference 
is larger for the exotic neutron star like material than for aluminum. The maximal difference is reached for the BS- 



Wc — Wf 
Wf 
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energy functional with roughly 3%. For increasing 5 the differences in three cases reach a maximum and afterwards 
decrease again. With increasing 6 the elastic energy in the configuration increases with respect to the gravitational 
rest mass energy. Thus, the more the energy of the configuration is dominated by the elastic energy the smaller is the 
influence of the choice of a relaxed state. In any case, what can be learned from Fig.|8]is, that for practical purposes 
one can safely assume that the metric of the relaxed state is flat. 



V. CONCLUSION 



We have discussed in this work the spherically symmetric body in relativistic elasticity for two different stored 
energy functionals. We find that the BS-fimctional corresponding to the classical Kirchhoff-St. Venant materials 
and the KM-functional have entirely different behavior for large deformations even though they agree for small 
deformations. The BS-functional gives rise to a mass-radius diagram which qualitatively is very similar to the one 
found by Karlovini and Samuelsson in | |T2|| . They obtain this diagram for a stiff ultra-rigid equation of state in 
the Carter-Quintana high-pressure formulation. They find that the decreasing branch is unstable. We can confirm 
this numerically and we even see indications of another region of configurations. This is an indication that the BS- 
functional gives rise to an increasingly stiff equation of state quite in contrast to the KM-fimctional for which the 
equation of state becomes increasingly soft. The result of this softness can be seen in qualitatively very different 
behavior of the mass-radius diagram which shows a spiral which approaches a limit point for large deformations. 

We looked at these functionals for two different materials, the 'every day' material aluminum and an artificial 
exotic material. While the sizes of the individual configurations are very different the qualitative behavior is very 
similar in both cases. 

In order to analyze in more detail the features in these configurations and in particular the stability properties of 
the different branches of the mass-radius diagrams it might be advantageous to formulate the problem not as an 
initial value problem as we have done here. Instead of specifying the central compression and integrating outwards 
to (possibly) find the boundary of a configuration one would instead set up a boundary value problem on the body 
subject to the boundary conditions imposed by the symmetry requirements in the center and the vanishing of the 
radial pressure on the boundary. Steps in this direction have already been made by Losert |9]. 

Of course, our considerations are to a certain point academic because any real material will break at already quite 
moderate deformations compared to the ones we have used. But we feel that such questions of principle may shed 
some light on the differences between the various possible choices and therefore on the justification of assumptions 
made when relativistic elasticity is used for real problems. 
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APPENDIX A: THE SYMMETRIC ENERGY-MOMENTUM TENSOR FOR THE KIJOWSKI-MAGLI ACTION 

The action for the Kijowski-Magli formulation of elasticity is | |2.15| | 

n(^ + ]^Aop- + ^-detgd^x. 

The energy-momentum tensor for this action is obtained (with our conventions) by variation of ^ with respect to 
the inverse metric 



r - o 



■ Pgab- (Al) 



The energy density p is specified in terms of the variables a = So" and j6 = jSa^'Sj," — ^{Sa")'^- The tensor Sj' was 
defined in ||ZT3t 

= log -t/)- (A2) 
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To find the variation in these variables we first need to compute the variation with g"^ in the tensor 

We sketch the calculations here without going too much into the details. We start with the variation of \2.5) to obtain 

Since, rp- = gah(-o''co^ we can now obtain the variation in n 

2— = Kb^g 

and with to" = nu" this yields the variation of and 

Su'' = ^ (^tkUaSg'^) M^ Sua = \ (ucUibg'^') Ua - gacU^Sg'^ . (A3) 

From this we find 

5{UaU^) = (ucU^Sg"^) UaU^ - gacU^uHg'^ = -hacU^ll^ 6g"^ (A4) 

and hence 

= -hacUdU^5g"^ - 'yacSg'K (A5) 

To compute 6x we use the formula 

TrlogA = logdetA (A6) 
which follows from the corresponding well-known equation 

det exp B = exp Tr B, 

valid for any matrix B. Thus, we have 

Set = -ijTrlogX = -i^logdetX = -i-^^^detK. 
2 ^ 2 * 2detX 

Using the formula 

^detK = detXTr(X-^M) 

we get the final result 

5cc = -^lr{K-^5K) = -^{K-^aSK^. (A7) 

The variation of /S foUows from the formula which can easily be derived 

^Tr/(A) = Tr(/'(A)M). 
We have j6 = ^ Tr(S^) — ^ (Tr S)^ so we first compute 

<5Tr(s2) = ijTr(log(X)2) = 1 Tr(log(X)X-iM) = -Tr(SX-i^X) = -Sa^K'^'SK,". 



With this result we obtain 



16^ = - IriSK'^SK) + ^TrSTr(K~^(5X) = - Tr{SK-^SK). (A8) 
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These formulas can be simplified as follows. Using the fact that Sa^'m" = and KJ'u" = we get 
Therefore, we have 

5^ = -\Ki,Sg'' 

and, similarly. 



Now we can write down the variation of the energy 



Sp-^-^p = n {AaScc + BSji) = -| [Aah.^ + BS,b] Sg''\ 
With this in hand we can now finally write down the energy-momentum 

Tah = p MflMfo + nOLAhal, + uBSab- 

APPENDIX B: ENERGY-MOMENTUM TENSORS 

Using the abbreviations x = F'e~^^~^'>] and y = F/r the energy-momentum tensor for the Beig-Schmidt theory 
is obtained from (|2.18t using the expression | |3.7| | for ej' in spherical symmetry. It has the following non- trivial 
components 

V = P = (po + f { (x' - 1)' + 2(/ - 1)2} + ^ { (x2 - 1) + 2(y2 - 1) }') , (Bl) 
Tji = -P = -i(xy2)x2 ((A + 2h)(x2 - 1) +2A(y2 _ l)) , (B2) 
Ti^ = = -Q = -J(x/)y' ({2\ + 2^){y^ - 1) + A(x2 - 1)) . (B3) 



The energy-momentum for the Kijowski-Magli formulation in spherical symmetry is obtained in a similar way 
from | |2.19l l once the tensor Sa^ and its invariants a and /5 are determined. Since S = — j log(l + e) and e is di- 
agonal in spherical symmetry this is straightforward: 

= -(logx)dr® a, - (logy) {dO (E}de + d(p(E}d^). (B4) 

Hence, we get 

a = -(logx + 21ogy), (B5) 

and 



1 



\2 , 0/1 n2 



1 



2 



/3 = - (logx)^ + 2(logy)^ --(logx + 21ogy)^ (B6) 



6 



With these expressions we obtain the non-trivial components of the energy-momentum tensor 

TflO = p = xy2 (po + ;i {(logx)2 + 2(logy)2} + ^ {logx + 21ogy}2) , (B7) 

Til = -P = -(xy2) ((A + 2y) logx + 2Alogy) , (B8) 

r2' = T33 = -Q= -(xy2)((2A + 2;i)logy + Alogx). (B9) 
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APPENDIX C: EXISTENCE AND REGULARITY OF THE SOLUTIONS 



Our goal here is to see whether the theorem by Rendall and Schmidt can be appHed to our situation. For easy 
reference, we cite the theorem here 

Theorem C.l (Rendall and Schmidt) Let V he a finite dimensional vector space, N : V ^ V a linear map all of whose 
eigenvalues have positive real parts, and G : V x (— e, e) — > V and g : (— e, e) ^ V smooth maps, where e > 0. Then, there 
exists 6 < e and unique hounded function u : {—S, 0) U (0, S) V which satisfies the equation 

du 

s — {s) + Nu{s)=sG{s,u{s))+g{s) (CI) 

Moreover, u extends to a smooth solution o/ (IClP on {—6,6). IfN, G and g depend smoothly on a parameter t and the eigenvalues 
ofN are distinct, then the solution depends smoothly on t. 



To analyze the behavior of the equations at the center we follow the paper by Park 112111 who has found that the 
equations for the Kijowski-Magli allow smooth regular solutions. We need the same result for the equations coming 
from the formulation of Beig and Schmidt. 

It is easy to see that a regular spherically symmetric scalar function u{r) on a spherically symmetric space-time 
can smoothly be extended as an even function of r to negative values of the radius. Thus, they depend smoothly on 
s = and their derivative at the center vanishes. Using s as the independent parameter instead of r our system (|3.19|l 
is 

2szv' + 3w = Anp, (C2) 
lsy' + y = xf (C3) 
IsP.x' = Py{y- fx)-sf^[w + AnP){P + p)-2{P-Q), (C4) 

where we have introduced / = (1 — 2sw)~^^'^. The solutions of this system will be regular only if x(0) = i/(0) = a, 
P(0) = Q(0) and w{Qi) = An/3p{0). Thus, we can write 

x = a + sxi{s), y = a + syi{s) 

and derive equations for the functions xi (s) and yi (s) . This yields the equations 

Isw' + 3w = 4:Tzp, 
f -I 

2sy[ + 3yi = a—^ h xif, 

2sxi+2xi = ^ (^i^fl+(i/i-/xi)) -£(ri; + 47rP)(P + p)-^^^. 

The desired form is 



2su' + Nu = sG{u,s) + g{s). 

The first equation is already in this form. In the second equation, the right hand side contains non-linear terms in w 
(via /) but these terms have no factor s in front. So we expand 

/ = (1 - 2swy'^^^ = i + sw + e{s^) 

and get 

2sy\ + 3i/i = aw + xi + ff{s). 

Thus, this equation is in the desired form. Note, that it is not important for this analysis to know the ^(s) exactly. It 
is enough to know that they are well behaved at the origin. The ^(s) implies that they have the desired factor of s in 
front. 

The third equation is more difficult.Here, the right hand side of the equations contains non-linear terms in w, Xi 
and j/i through the dependence of / and the functions P, p and Q. Therefore, we have to expand the right hand side 
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around s = and see whether we can isolate constant linear terms and non-linear terms proportional to s. So we 
write 

P = Pix,y) = P{a + sxi,a + syi) = P{a,a) + s{xiPx{a,a) +yiPy{a,a)) + &{s^) 

and similarly for all the other functions on the right hand side. We write P{a,a) = P^ and similarly for the other 
functions. First, we observe that we need to have 

P° - = and P° ^ 

for the right hand side to be regular at all. Expanding the right hand side as indicated we get 

This is the desired form and we can collect the coefficients of the linear terms in the matrix N 



N 



/ 3 

-a 3 

V — ^-w^ 




For the theorem to apply the eigenvalues of this matrix must be positive. This leads to the following conditions: if 
P° > then also 

6P0 + 2P0-3QS-Q0 >0, (C5) 
7p0 + pO _ 2q0 > 0, (C6) 

while for P^ < the above expressions must also be negative. 

Now we have reduced the question of regularity at the center to a condition for the energy functional which 
specifies the elastic properties of the material. It is straightforward to check that these conditions are satisfied for 
both the Beig-Schmidt and Kijowski-Magli stored energy functionals. 
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